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ABSTRACT 

We calculate light-curves for Gamma-Ray Burst afterglows when material ejected in the explosion is confined 
to a jet which propagates in a medium with a power-law density profile. The observed light-curve decay steepens 
by a factor of T 2 when an observer sees the edge of the jet. In a uniform density medium the increase in the 
power-law index {(J) of the light-curve as a result of this edge effect is ~ 0.7 and is completed over one decade 
in observer time. For a pre-ejected stellar wind (p oc r~ 2 ) (3 increases by ~ 0.4 over two decades in time due to 
the edge effect and the steepening of the light-curve due to the jet sideways expansion takes about four decades in 
time. Therefore, a break in the light-curve for a jet in a wind model is unlikely to be detected even for very narrow 
jets of opening angle of a few degrees or less, in which case the lateral expansion occurs at early times when the 
afterglow is bright. 

The light-curve for the afterglow of GRB 990510, for which an increase in f3 of approximately 1.35 was ob- 
served on a time scale of 3 days, cannot be explained only by the sideways expansion and the edge effects in a 
jet in a uniform ISM - the increase in j3 is too large and too rapid. However, the passage of the cooling or syn- 
chrotron peak frequencies through the observing band at about 0.1-1 day together with jet edge effect explains 
the observed data. The jet opening angle is found to be ~ 5° and the energy in the explosion to be less than about 
10 50 erg. 



Subject headings: gamma rays: bursts - gamma-rays: theory 

1. INTRODUCTION 

In a recent paper Chevalier and Li (1999) pointed out that 
some of the GRB afterglow light-curves are best modeled when 
the density of the circum-burst medium is taken to fall off as 



r (this is referred to as the wind model). These afterglows 
show no evidence for a jet, i.e. their light-curves follow a 
power-law decline without any break. This is puzzling since 
collimated outflows are expected in the collapsar model for 
GRBs (MacFadyen, Woosley & Heger 2000). We offer a possi- 
ble explanation for this puzzle by showing that the light-curve 
resulting from the interaction of a jet with a pre-ejected wind 
falls off as a power-law whose index changes very slowly with 
time. 

We carry out a detailed modeling of the multi-wavelength 
afterglow flux data for GRB 990510, which provides the best 
evidence for a jet propagation in a uniform density medium 
(Harrison et al. 1999, Stanek et al. 1999), to show that ef- 
fects associated with a finite jet opening-angle are insufficient 
to explain the observed rapid steepening of the light-curve. 

In §|2| we calculate the propagation of a jet in a stratified 
medium and in we describe the calculation of the syn- 
chrotron emission and afterglow light-curve. 

2. DYNAMICS OF EXPANDING JETS 

The dynamical evolution of jets and its synchrotron emission 
have been previously investigated by a number of people, e.g. 
Rhoads (1999), Panaitescu & Meszaros (1999), Sari, Piran & 
Halpern (1999), Moderski, Sikora, & Bulik (2000), Huang et 
al. (2000). The evolution of the Lorentz factor (T) can be cal- 
culated from the following set of equations 
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where 8 is the half-opening angle of the jet, Mo, and Tq are the 
initial mass and Lorentz factor of the ejecta, M\ is the swept-up 
mass, p(r) — Ar~ s is the density of the circum-stellar medium, 
and f — c/c s is the ratio of the speed of light to that of the jet 
sideways expansion; / is a parameter of order unity whose ef- 
fect can be absorbed in Tq, and which has little effect on the 
light-curve. O is the angle between the velocity vector at the 
jet edge and the jet axis (in the lab frame) and is determined 
by the modification of particle trajectory due to the sideways 
expansion. The last equation above expresses the conservation 
of energy and it applies to an adiabatic shock when the heating 
of the original baryonic material of rest mass Mq by the reverse 
shock is ignored. 

The above equations can be combined and rewritten in the 
following non-dimensional form which is applicable for rela- 
tivistic as well as non-relativistic jet dynamics 
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is the deacceleration radius. E = MqTq is the energy in the 
explosion and 9 is the initial half-opening angle of the jet. The 
above equations show that for the wind and the uniform ISM 

]/4 3/8 1 

models T oc t ob J & t ob J , respectively, as long as T ^ 9q , 
where t a i, s = J dt(l — v) is the observer time, t being the lab 
frame time and v the jet velocity in units of c. 

Equations (Q) and (g) are solved, subject to the boundary con- 
ditions 2/1=2/2 = 1 for x <C 1, to determine T and 9 as func- 
tions of r. For a relativistic jet with Q = 9, i.e. fluid velocity in 
the radial direction, equation (|5|) reduces to 
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The solution of the equations (Q) and (|j) is a two-parameter 
family of functions, however in the relativistic case the solution 
depends only on the product 9qTq. 

One can solve equation (0) approximately, ignoring the very 
early time behavior, to determine the time when the sideways 
expansion alters significantly the jet dynamics. The two rela- 
tions in equation (Q) can be combined to yield a first order dif- 
ferential equation for 2/12/2 = y which is given by 
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with T) = /(3 — s)(9qTo), a constant, and £ 
An approximate solution to this equation is 
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Thus, y oc T9 decreases monotonically with radius or time. 
The transition to jet sideways expansion starts when the two 
terms in the above equation become equal, i.e. £ ~ (r//16) 2 , 
and lasts for an interval in £ for which y decreases by a fac- 
tor of ~ 3, or x increases by a factor of ~ 3 3 ^ 3 ~ s \ The 
Lorentz factor continues to fall during the transition by a fac- 
tor of a few. Therefore, the transition time divided by the time 
at the start of the transition (in observer frame), during which 
ai = — rfln(T — l)/d\nt {, s increases from (3 — s)/(8 — 2s) to 
approximately 1/2, is approximately 9x3 3 ^ 3_s ^. The solution 
to 2/1 and 2/2 can be obtained by inserting the expression for y 
into equation (0). However, 2/1 and 2/2 determined this way have 
much larger error than y and should not be used for any serious 
calculation. 

We solve equations (0) and numerically and show the re- 
sults for x(t i, s ) and a\[t i, s ) in Figure 1. Note that the change 
to a\ from one asymptotic value, corresponding to spherical 
shell expansion, to another, when sideways expansion is well 
underway, takes a long time; the ratio of the final to the ini- 
tial time for a change in ot\ of 0.1 for a uniform ISM is ~ 10 2 
whereas for s = 2 the ratio is 10 3 . For the parameters chosen 
here ol\ — 0.5 when T is of order a few. In the non-relativistic 
phase of the jet expansion ai\ — 1.2, as for a Sedov-Taylor 
spherical shock wave. 

3. SYNCHROTRON EMISSION FROM RELATIVISTIC JETS 

The synchrotron spectrum in the co-moving frame is taken to 
be a sequence of power-laws with breaks at the self-absorption, 



synchrotron peak, and cooling frequencies, as presented in Sari, 
Narayan & Piran (1998); these frequencies can be found in eg. 
Panaitescu & Kumar (2000). All of our numerical results, un- 
less otherwise stated, are obtained by integrating emission over 
equal arrival time surface. Ignoring the radial structure of the 
jet, the flux received by an observer located on the jet axis is 
given by 
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where P' v , is the co-moving power per frequency at 
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We ignore the angular integration when discussing the ana- 
lytical calculation of the observed flux and its power-law de- 
cline with time. The observed flux at a frequency that is greater 
than both the cooling frequency, v c , and the synchrotron peak, 
v m , is proportional to 
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At early times when T 9 1 and r oc t 

decays as t~ b ^ p ~ 2 ^ 4 . At late times when T9 < 1 the power- 
law index for the flux (3 = — d In f v /d In t b s — (4 — s) [ai (p + 
2) — l]/2 + sp/4 + 0:2, where 0:2 = — 2d\ny /d\nt b s - 

There are two effects that determine the evolution of j3. One 
of them, the edge effect, is purely geometrical and results from 
the angular opening ~ of the relativistic observing cone 
becoming larger than the jet opening angle 9, i.e. the observer 
"sees" the edge of the jet. The increase to (3 resulting from it 
is oii 5; (3 — s)/(4 — s); a 2 decreases with time and therefore 
the jump in (3 is smaller for larger #o- The dimensionless time 
for (3 to increase by 012 depends on the angular position of the 
observer w.r.t. the jet axis and is approximately the ratio of the 
time when the observer sees the far edge of the jet to the time 
when the near side of the jet becomes visible. This time is given 
by 

( 8 -aa)/(3-.) r 1 + p i/2i( 8 - 2 *)/( 3 -s) 

Rt 
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where P^q is the probability that the observer lies within an an- 
gle (j) of the jet axis. For P 0O = 0.25, R tc is 18.7 (81) for 
s = (2), and during this time (3 increases by approximately 
0.7 (0.4). The dependence of R tc on 0o becomes much weaker 
when the emission is integrated over equal arrival time surface 
(Figure 2). This is because the effect of angular integration is 
to smear the jet-edge by an angle 1/r ~ 6q/2, which sets the 
minimum value of R ts to be about 10 (10 2 ) for uniform (wind) 
models. 

The other effect which leads to a steepening of the afterglow 
decay is dynamical and is caused by the lateral spreading of 
the jet. During the relativistic phase the increase to (3 from the 
sideways expansion is 5(3 = (p + 2)(4 — s)8a\/2 + 5a.2\ 
Sa± and 8ct2 can be read from Figure 1. Since ct\ does not 
asymptote to 0.5 so (3 ^ p during the relativistic sideways ex- 
pansion of the jet. 1 The value of (3 does, however, approach 

1 It should be noted that the asymptotic behavior /3 — > p for s = (Rhoads 1999) is achieved only for extremely narrow jets (80 < 1°). so that the jet remains 
relativistic for a sufficiently long time after it starts expanding sideways. It nevertheless serves as a useful, quick way of estimating p approximately from the late time 
light-curve, when /3 is no longer increasing. 
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p because Sat « 1/(8 — 2s) sometime before the jet becomes 
non-relativistic and a 2 ~ at this time, thereby giving j3 « p 
(see eq. [[j~2|] and Figures 1 & 2); /? can exceed p, as can be seen 
in Figure 2, however the decrease in a 2 during the mildly rela- 
tivistic phase prevents (3 from getting much larger than p. This 
result can be extended to any observing frequency v > v m after 
an appropriate modification of equation < fl2[ > . For instance, to 
consider the case of v c > v > v m the right side of the equation 
should be multiplied by a factor of (t b s r 2 )( 1_3s / 4 ), which has 
little effect on the evolution of (3. The time scale for the in- 
crease in (3 due to sideways expansion is of order 10 2 (10 3 ) for 
s=0 (2) (see Figure 2). Therefore this effect is smaller than that 
resulting from seeing the jet edge, and it extends over a much 
longer time. 

To conclude, we wish to emphasize that for most jets prop- 
agating in a uniform ISM we are likely to see an increase to (3 
of only 0.6-0.9; the remainder of the increase takes place on a 
long time scale, and thus is hard to detect. For jets in a windy 
medium, s = 2, j3 changes by less than about 0.5 and the tran- 
sition time R te ~ 10 3 . Such a gradual increase to the afterglow 
light-curve power-law index is extremely difficult to detect (see 
Figure 2). For instance, if the edge of the jet becomes visible at 
tabs ~ 1 day, the difference in the optical flux at the end of 10 
days with and without jet is ~ 0.25 mag, which can be easily 
missed. Thus, the GRBs studied by Chevalier and Li (1999), 
which show evidence for the wind model, could in fact have 
had a collimated ejection of material. 

3.1. The Afterglow of GRB 990510 

The optical emission of the afterglow of GRB 990510 was 
measured in the V, R and I bands between 0. 15 and 7 days after 
the burst and showed the power-law index of the light-curve, (3, 
to have increased from 0.82 ± 0.02 to 2.18 ± 0.05 (Harrison et 
al. 1999)orfrom0.76±0.01to2.40±0.02(Staneketal. 1999) 
during a dimensionless time i? tc « 30 which, as described pre- 
viously, is not possible to obtain through the effects of the jet 
sideways expansion alone. Therefore there must be some con- 
tribution to the light-curve steepening due to the passage of one 
(or both) of the spectral breaks: the synchrotron peak v m and 
the cooling frequency v c . 

In Figure 3 we show a comparison between the light-curves 
of GRB 990510 in the V, R, I bands and the 8.7 GHz radio 
data, with a model where the cooling frequency v c crosses the 
optical band at t Q b s ~ 1 day. The steepening of the light-curve 
has little dependence on the observing band because the ratio of 
the largest to the smallest optical wavelength is ~ 1.5. More- 
over, the integration over angle spreads in time the steepening 
of (3, making it nearly achromatic. An increase of (3 by ~ 0.8 
is caused by the jet edge and the sideways expansion, and an 
increase of 0.25 results from the passage of v c through the ob- 
serving band. A further increase of (3 of ~ 0.15 is caused by 
the passage of v m through the observing band at t a i, s ~ 0.03 
day (see lower panel of Figure 3); the transition time for [3 to in- 
crease by ~ (3p— 1)/4 due to the v m crossing is about a decade 
in the observer frame as a result of integration over equal arrival 
time surface, hence one should be careful in deducing p from (3 
at early times. All these together give rise to a light-curve that is 
consistent with the data. The model is also consistent with the 



HST V^-band observation carried out at about a month after the 
burst (Fruchter et al 1999). The parameters for the fit are given 
in the caption for fig. 3 which yields the energy in the burst to 
be 2xl0 49 erg. Correcting for the radiative losses the energy 
in the burst increases by a factor of a few to < 10 50 erg. We 
estimate the uncertainty in model parameters by varying them 
in such a way that the numerically calculated light-curve lies 
within 3-<7 of the observed data points. We find the uncertainty 
in the jet opening angle and the burst energy to be a factor of 
two, and e e , n and eb are found to be uncertain by factors of 
about 4, 40 and 7 respectively; we note that the radio observa- 
tions are very important in constraining the model parameters. 
The electron index p is constrained by the observed (3 before 
and after the ~ 1 day break; the error in p is ~ 5%. 

The optical emission of the afterglow of GRB 990510 can 
also be explained by a model where the synchrotron peak fre- 
quency crosses the observed band at ~ 0.1 day. Its effect on 
f3 persists for up to ~ 1 day and yields an increase of (3 of 
~ 0.5 during the early observations. The parameters for the 
second model differ from the one described above (see Figure 
3) somewhat. In particular, e e is larger by a factor of two, the 
energy per solid angle is smaller by a factor of two, and 9q is 
- 20% larger. 

4. CONCLUSIONS 

One of the main results of this work is to show that afterglows 
from well collimated Gamma-Ray Burst remnants going off in 
a medium with density decreasing as r -2 show little evidence 
for light-curve steepening due to jet edge and sideways expan- 
sion. This could explain the lack of breaks in the afterglows of 
GRB 980326 and GRB 980519, which Chevalier & Li (1999) 
found to offer support for the wind model. Jets can perhaps be 
detected by the measurement of time dependent polarization. 

In a collimated outflow the sharpest break in the light-curve 
is produced in a uniform density circum-stellar medium, and is 
associated with the edge of the jet coming within the relativistic 
beaming cone (the edge effect). The magnitude of this break is 
~ 0.7 (0.4) for a uniform ISM (wind model) and occurs over 
about 1 decade (2 decades) in time. Further steepening of the 
light-curve, associated with the sideways expansion of the jet, 
occurs on a much longer time scale of R te ~10 2 (10 4 ), i.e. 
weeks to months. 

The power-law index for the light-curve of GRB 990510 in- 
creased between days 0.8 and 3 by about 1.35. This is too large 
and too fast to result from jet edge & sideways expansion ef- 
fects. However, the observations can be explained if either the 
cooling or the synchrotron peak frequency passed through the 
observing band at about 1 or 0. 1 day, respectively. Models that 
are consistent both the optical and radio data of this afterglow 
have an opening angle of ~ 5° and energy in the explosion is 
< 10 50 erg (see Figure 3). 

For the afterglow of GRB 990123 the power-law index of the 
light-curve increased by 0.55 between days 1.5 and 3, which 
can be explained by the edge effect alone (Meszaros & Rees 
1999). 

We thank Peter Meszaros and Vahe Petrosian for useful dis- 
cussions. 
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FIG. 1. — The main panel shows a\ = — dln(r — 1) / d\nt ^ s for uniform ISM s = (thin continuous curve) and the wind model s = 2 (thick continuous 
curve), and a?2 = —2dln(8F) /dlnt i, a for s=0 (thin dashed curve) and s=2 (thick dashed curve). The inset shows the evolution of the jet radius. The symbol 
marks the time when 6T = 1, and the symbols 5 (2) denote the time when 7 = 5 (2). For this calculation we took the energy per unit solid angle to be 3 X 10 53 erg 
sr" 1 and 6 = 1/30 rad. The density of the uniform ISM is 1 cm" 3 , A = 5 X 10 11 gem" 1 for s = 2, and T = 300 (i.e. r 6» = 10). 
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FIG. 2. — The panel shows /3 = — din /„/ In t j, s for v > max(^ m , v c ) and </>o = 9o/V5. The thin & thick dashed lines are for s = & 2 respectively; the 
observed flux in these cases was calculated without proper angular integration over the jet surface. The sharp increase to the value of /3 seen in the dashed curves 
arises from the edge effect described in the text. The thin and the thick continuous curves are for s = & 2 respectively and these calculations included integration 
over equal arrival time surface; note that this smoothes out sharp changes in f3. For all of these calculations we took the energy per unit solid angle to be 3 X 10 53 
erg sr — 1 , 9q = 1/30 rad (BoTo = 10), p = 2.5, the density of the unifoon ISM is 1.0 cm -3 and A = 5 X 10 11 gem -1 for s = 2. The results shown here are 
independent on the sideways expansion speed (i.e. the parameter / in eq. [fc||), as long as this speed is relativistic. 
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FIG. 3. — Comparison between the observed and the theoretically calculated light-curves for the afterglow of GRB 990510 in the radio and the optical hands Most 
of the data is taken from Harrison et al. (1999); the optical data was supplemented with observations reported in GCN Circulars (http://gcn.gsfc.nasa.gov/gcn/gcn3) 
by Beuermann et al. 1999, Marconi et al. 1999, Pietrzynski & Udalski 1999; the latest measurements in the V-band were taken with the HST (Fruchter et al. 1999). 
The model light-curves are calculated for an observer located on the jet axis, which gives the fastest decline of the light-curve. The jet has energy per solid angle 
_B/(7r6»g) = 1.2 X 10 52 /4tt erg sr and a half-angle do = 0.085 rad. The electrons acquire e e = 0.3 of the internal energy after shock acceleration, the magnetic 
field energy is eb = 0.04 of that of the shocked gas, and the electron index is p = 2.2 . The external medium is homogeneous with n = 0.23 cm~ 3 . The redshift 
of the source is z = 1.62. The cooling frequency passes through the observing window at t (, s = 1.2 day steepening the afterglow light-curve while the sideways 
expansion is effective. The inset of the upper panel shows the 8.7 GHz emission. The lower panel shows a comparison of the numerically computed power-law index 
(/3) for the decline of the afterglow of GRB 9905 10 and the observed one, as obtained by the fitting formula given in Harrison et al. (1999). 



